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RADIOMETER DATA PROCESSING 
IN THE HAYSTACK ANTENNA POINTING SYSTEM 


L. INTRODUCTION 

Since Radiometric methods are used to boresight the Haystack System, the 
Univac 490 antenna pointing system includes a radiometer data processing program 
which operates in real time concurrently with the antenna pointing programs. 

The program implements a technical approach which is matched to the radio- 
metric measuring equipment* provided by Messrs. Meeks and Weinreb of Group 31. 
This approach was jointly formulated. 

The main input to the program is comprised of 30-bit data words passed from 
the radiometric equipment to the computer via channel 5. Each word contains 
identification bits allowing the computer program to distinguish between two different 
types of data: “auxiliary” or radiometric’. The processing of the auxiliary data 
consists of a simple check to determine if data associated with a particular identifi- 
cation falls within specified bounds for that identification. This data is then logged 
on the high-speed printer, along with its associated upper and lower limits, and an 
indication if the data is outside of the specified limits. Processing of the radiometer 
data is more complex and is described in the following sections; in general, the 
primary Output is a running measure of radiometric source temperature. The results 
of this processing are both logged and plotted on the high-speed printer in real time. 
lk. INPUTS 

A. Via Channel 5. The bit layout of the incoming words is shown in Figs. 
land 2. The incoming word is either auxiliary data (identified 1 through 50) or 
radiometric data (identified 51 through 53). ” The data entering the system can be 
thought of as existing in blocks. A block of auxiliary data (A data) can be up to fifty 


words. The auxiliary data train may be somewhat shorter in length and may consist 


* ‘The radiometric Equipment is described in TR-365, H. G. Weiss, 15 September 
1964. 


* Note that an input word cannot contain 0 (30-bits of 0) or, —0 (30 bits of 1) by 
virtue of the identification field. 


of any arbitrarily selected subset of the fifty words. the only restriction being that 
lower numbered [. D.'s must be sent first. A block of radiometric data (R data) is 
three words long, the first being a status word. (This can be readily modified to 
accommodate blocks of R data up to fifty words long.) The bit layout of the status 
word is shown in Fig. 2. The status word indicates one of four possible statuses - 
Base, Calibrate, Observe, or Stop. The next two words of the radiometer data 
block contains the data from receivers l and 2. [Uf the incoming word starts a block 
of auxiliary data or radiometer data, the right ascension and declination of the point 
on the celestial sphere corresponding to the direction of the antenna (saved during 
last integral second) is stored along with sufficient data to determine time to the 
nearest 4ms. A sense switch (SS3) is included in the radiometer equipment. When 
set by the experimenter, SS3 results in the printing of an asterisk at the beginning of 
an R data output line. 

The data rate is selectable from the radiometer equipment. The antenna point- 
ing system is interrupted upon the arrival of a 30 bit word. Incoming data is stored 
in a circular buffer consisting of 125 registers which can accommodate up to 20 
complete R data blocks. This buffer is cleared once each frame (nominally 2 sec. ). 
The data flow to the computer is asynchronous with any other antenna pointing system 
function. 

B. Keyboard Inputs 

When the radiometer program is initialized, the experimenter Is given the 
option to change any or all auxiliary data limits and any or all calibration constants 
(those used in equations 3 and 4). The auxiliary data limits are entered as decimal 
integers (— 9999 through +9999). The calibration constants are entered as decimal 
numbers. 

The experimenter may also elect to intersperse “comments” with his program 
output. Comments are typed in via the console and can be up to 80 characters per 
line, terminated by a carriage return. Figure 8 shows a sample keyboard input. 

Ge Common Storage Input 


The following common storage registers are used. 


TRUE TIME = gmt 


CAZIM > corrected azimuth 
CELEV “ corrected elevation 
ASTRORA : displayed right ascension 
ASTRODEC displayed declination 


Il. DATA PROCESSING 

The processing of data is discussed under the headings of Radiometric Data 
and Auxiliary Data. 

A. Radiometric Data 

There are three distinct categories of radiometric data, Calibrate, Base and 
Observe. A Calibrate run can be started any time, may be of arbitrary length as 
chosen by the experimenter, and is terminated upon the arrival of data marked as 
Base. A Base run can be started only after a Calibrate run, is again of arbitrary 
length, and is terminated by incoming Observe Data. Each base or calibrate run is 
considered as a single integration period, whereas the observe run is broken up into 
a multiplicity of integration periods; the length of the Observe integration period (in 
Blocks) is specified by the status word. Over each integration period of calibration, 
base or observe data, the following quantities are determined; it should be noted that 
a particular quantity will be produced for each of the two radiometers. 
The average signal is 


l N 
Ney (1) 


where N is the number of blocks and r is the data point. 

The estimated standard deviation is 

N 2 

yw. £X, 2 (2) 


ul i R 
V dadsoy - NoT fox N > 2 


Subsequent processing will depend on the particular radiometric data category. 


.. 


l. Calibrate 
For each calibrate period, the following quantities are defined for later proces- 


Sing and on-line printing. 


B is the average signal during the Calibration period. 
So is the estimated standard deviation. 


Nis the number of blocks in the Calibration period. 
Cc 


pe Base 


For each Base period, the average signal and standard deviation define the 
following quantities. 


Ry is the average signal during the Base period. 


Sb is the estimated standard deviation. 


Ny is the number of blocks in the Base period. 


The following quantities are then computed for each receiver. 
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where X and Y are constants called TCAL and TBASE. 
Bs Observe 

It has been noted that the integration period used for the Observe category is 
specified in the status word. In the event of a change in the number of blocks per 
integration period, data accumulated up to the time of change will be included as the 
first samples of the new interval. [nthe event of a change in the status word to a new 
data category, the accumulated data will be dropped. 

For each Observe interval, the following quantities are defined. 

R is the average signal during the integration period. 

55 is the estimated standard deviation. 


2a is the number of blocks in the integration period. 


The following quantities are computed for each receiver 


= * R —R (R, — Ry) 
4 b 
_ 4 
AT = y—y S (4) 
Cc b 
4, Stop 


The stop category is used to make the program "idle". If the previous cate- 
gories were either CAL, BASE or STOP the program rejects this block and proceeds 
to look for more data. If the block immediately preceding this stop block was of the 
OBSERVE category, the accumulation registers are cleared (this discards the last 
incomplete integration period) and spaces the line printer one line. The program 


then proceeds to look for more input. 


B. Auxiliary Data 


For each Aj, an upper and lower limit is stored. These limits may be changed 
individually via the keyboard during the initialization phase of the program. 

Processing of the data is a simple check to determine if the data for an A. 
falls within the specified limits. 
IV. PROGRAM OUTPUTS 


A. Real Time Outputs 


The real time output of data is via the high-speed printer. Figure 3 shows a 
sample of data output on the printer after a Base run. ‘This output is triggered at the 
end of a Calibration-Base sequence and is printed at the top of a page. The duration 
is in blocks and is mt. and Ny whereas DELTA CAL (1), DELTA CAL (2), DELTA 
BASE (1), DELTA BASE (2), TEMP (1) and TEMP (2) are found from Eq. (3). The 
calibration constants are those used in all computation. The values for azimuth, 
elevation, right ascension and declination are those for the antenna position at the end 
of the BASE category run. 

Figure 4 shows the format for the logging and plotting of the OBSERVE data. 
The column headings and scaling information (showing the range of the plot) are 
printed as the second line on each page. Time, right ascension and declination are 


those for the midpoint of the integration period. In the event the scale is changed, 


the plot shows a discontinuity allowing the new scale to be printed. The scale is 
selected via the status word. A change in scale (and hence range) does not interrupt 
the processing of data. The symbol for receiver 1 is "X" and for receiver 2 is 
"O". When both receivers have the same value, a single symbol "1" is plotted. If 
a quantity exceeds the plotting range, it will be plotted as the closest value within 
the selected range. A blank line signifies discarded data. 

Figure 5 shows the format for Auxiliary Data. The column headed by E is 
used to denote that the value for a particular A; lies outside its limits. The asterisk 
character will denote this phenomenon. The auxiliary data printout will always be on 
a separate page. 

Be Emergency messages appear on the high speed printer. The printer is 
spaced so that the message can be read without using the line feed. There are two 
emergency messages, one calling for a calibrate sequence, the other asking fora 
Base run. A calibrate sequence is requested if the experiment is not started with 
Calibrate data. A Base run is called for if an experimenter tries to follow a Cali- 
brate run with Observe data. 

V. MAGNETIC TAPE RECORDING 


In keeping within the antenna pointing system recording philosophy the radio- 
meter program will prepare the data records and the system recording program will 
perform the actual tape writing (binary records) and error checking. These records 
will then appear on the system recording tape along with any other system recordings. 
The radiometer program prepares two types of records, data and comments, each 
having a unique [D. 


A. Comments 


A comments record will always be 18 words (30 bit) long. The first word will 
contain the field data coded characters "RDMTI". The second word will indicate a 
writing parity for the previous record. (0 = noerror). The next 16 words will 
contain the tield data coded characters that were typed in as comments. One of 
these records will be written each time the experimenter terminates a comment line 


via the carriage return. 


B. Data Records 


The Data records will always be 152 words long. The first word contains the 
field data coded characters "RDMTR". As in the case of the comments record, the 
second word is used by the recording program and sets a non zero value (30 bits) if 
a writing error occurred for the previous record. At present 4 types of informa- 
tion are recorded, each having their own sentinel. The four types of information may 
be intermixed on the recording in any order. If for any type, the 152 word limit is 
reached the rest of the words will appear at the beginning of the next radiometric data 
record. 

‘8 Radiometric data sentinel is 77777 00001 

A sentinel of 77777 0000! indicates the next seven words are from a radio- 

metric data block. 
word is the output azimuth buffer control 
word is true time in days witha B of 27 


is radar elevation from the encoders. B of 19 


] 
2 

word 3 is radar azimuth from the encoders. B of 19 
word 4 
5 


word is the status word (see Fig. 2) 
word 6 is from receiver 1 (see Fig. 1) 
word 7 is from receiver 2 (see Fig. 1) 
The first four words are recorded upon the arrival of a status word. Words 
1 and 2 permit a calculation of time to the nearest 4ms. as follows: 
word 2 gives time at the beginning of current computer frame. 
The upper half of word 1 gives the final location for azimuth 
output for this frame, while the lower half gives the next loca- 
Gon. thareiere final address - current address gives the 
fraction of the frame that has elapsed since truetime. 
is Auxiliary data sentinel is 77777 00002 
The first four words following this sentinel are the same as for radiometric 
data. They are recorded when the first word of an Auxiliary data scan is sensed. 


These words are followed by as many (up to 50) auxiliary data words as are sent to 


the computer. 


3. Calibration sentinel is 77777 QOO003 
Four calibration constants follow this sentinel: 


word | is Y ,(TBASE )) BIS 
word 2 is Y.(TBASE.,) B15 
word 3 is X (TCAL )) B20 
word 4 is X,(TCAL,) B20 


This data is recorded during the initialization phase of the program at the 


beginning of an experiment and any time the constants are changed. 


4. Right ascension and declination and range are recorded after 
the sentinel 77777 QO004. The angles are recorded in revolutions with a B of 27. 
Range is recorded in earth radii witha B of 22 if positive. If range is negative, 
the compliment of the range in astronomical units with of B of 22 is stored. If 


range is O an infinite range has been assumed. 


VI. PROGRAM DETAILS 
A logical flow diagram of the computer program is shown in Figs. 6 and 7. 
The three main sections of the program are initialization, interrupt and working. 
A. Initialization 


The initialization section performs the following function: 


L. Clears the accumulation registers (Nr), etc. ). 

5a Clears input buffers and sets control for selecting next word in 
buffer. 

a Sets control for storing incoming data into work blocks. 

4, Clears output line count. 

3. Gives experimenter the option to change A data limits. 

6. Types on console the calibration constants and allows the exper- 


iments to change these constants. 

Fe Gives the experimenter the opportunity to type comments and have 
them appear on the high-speed printer, and on the magnetic tape 
recording. (Although this portion of the radiometer program is 
in the initialization section, the facility for typing comments 


ren.ains active throughout the experiment. ) 


The initialization of the radiometer program is accomplished through the con- 
sole typewriter. Figure 8 shows a page copy of a sample operator/program communi- 
cation. In the example the experimenter changes a few A data limits and the calibration 
constants. 


B. Interrupt 


The antenna pointing system is interrupted upon the arrival of a 30 bit word 
from the radiometer equipment. Control is then passed to the interrupt portion of the 
radiometer program. ‘The interrupt section first saves the operational registers 
then stores (in BUF) the incoming word for later processing. If the incoming word 
signifies the start of a new block a 30-bit minus 0 is stored. Time information is 
then stored. The operational registers are restored, channel 5 is reset to provide an 
interrupt upon the arrival of the next word and control is returned to the antenna point- 
ing system. 

The beginning of a block is determined from the identification of the word and a 
control indicator called LASTADIND. The beginning of a Radiometric data block is 
merely a word with an [ID of 5l. The beginning of an Auxiliary data block is sensed 
if the [Dis less than 51 and if the LASTADIND is set to a non zero value. If the 
beginning of an A data block is sensed, LASTADIND is set to zero. This indicator is 
reset by the working section of the program. 


oF Working Section 


The working section operates as a data processing program immediately follow- 
ing the azimuth buffer chain. A routine called "get next block” is used to process the 
buffered data. This routine checks the circular buffer for any data (non-zero words) 
and sets up a work block which may be 52 registers long. The first word of the 
work block dictates the setting of the work block indication (0 = R/data, 1 = 
A/data. Inthe event of R data the plotting scale is sensed - if a new scale is 
indicated, it is printed and the line count is increased by one. Control is passed to 
the working section if a complete block is found - if not control is passed to the 
master control program. (See Fig. 7) 

The processing of an A data block results in a printout as shown in Fig. 5. It 
is noted that if a particular A(1l) is blank no data was received for that ID. The 
LASTADIND is then set to accept the next A data scan and control is passed to the 


master control program. 


In the event the work block is R data, the data status is sensed. ‘The initial 
data status is stop. Figure 7 shows the switching logic which controls further proces- 
sing of the data. Note that the first R data that can be accepted must be of the 
CALIBRATE mode and that this must be followed by a base run. When a block is 


accepted the data are included in the following calculation. 


- 29 

~T. BO max value 2 

Er BO double register - 
8) 

N BO max value 2 


When a particular mode is terminated, Eq. (1) and (2) are used to calculate. 


Bak it BLS 


* 
a B13 


The values for V, AV, AC as given by Eq. (3) are calculated upon the completion 
of a base run. They are also scaledtoa B of 15. 

Finally T and AT are calculated (witha B of 15) for each observe inte- 
gration period. These are also scaledtoa B of 15. The results are then recorded 


on the high speed printer via the system logging program. 


*It should be remembered that these values are calculated for each of the two 
receivers. 
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